Quadratic solitons as nonlocal solitons 
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We show that quadratic solitons are equivalent to solitons of a nonlocal Kerr medium. This 
provides new physical insight into the properties of quadratic solitons, often believed to be equivalent 
to solitons of an effective saturable Kerr medium. The nonlocal analogy also allows for novel 
analytical solutions and the prediction of novel bound states of quadratic solitons. 
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Quadratic nonlinear (or x^ 2%> ) materials have a strong 
and fast electronic nonlinearity, which makes them ex- 
cellent materials for the study of nonlinear effects, such 
as solitons 0. The main properties of quadratic soli- 
tons are well-known @ and both 1+1 3] and 2+1 
dimensional bright spatial solitons have been observed 
experimentally. Unlike conventional solitons, which form 
due to a self-induced refractive index change, the forma- 
tion of quadratic solitons does not involve any change 
of the refractive index. Thus the underlying physics of 
quadratic solitons is often obscured by the mathemat- 
ical model. Only recently Assanto and Stegeman used 
the cascading phase shift and parametric gain to give an 
intuitive interpretation of effects, such as self-focusing, 
defocusing, and soliton formation in x^ materials ||. 

Nevertheless certain features of quadratic solitons, 
such as soliton interaction, are still without a physical in- 
terpretation. Here we use the analogy between paramet- 
ric interaction and nonlocality and present a physically 
intuitive nonlocal theory, which is exact in predicting the 
profiles of stationary quadratic solitons and which pro- 
vides a simple physical explanation for their properties 
including formation of bound states. The nonlocal anal- 
ogy was applied recently by Shadrivov and Zharov to find 
approximate bright quadratic soliton solutions, but the 
nonlocal concept was not fully exploited to give a broad 
physical picture in the whole regime of excistence 0- 

We consider a fundamental wave (FW) and its second 
harmonic (SH) propagating along the z-direction in a 
quadratic nonlinear medium under conditions for type 
I phase-matching. The normalized dynamical equations 
for the slowly varying envelopes Ei_ 2 (x, z) are then Q 

id z Ex + d x d 2 x E x + E\E 2 exp(-i/3z) = 0, (1) 
id z E 2 + d 2 d 2 E 2 + E 2 exp(i/3z) = 0. (2) 



In the spatial domain d\K,2d 2 , dj>0, and the coordinate 
x represents a transverse spatial direction. The term 
d 2 Ej then represents beam diffraction. In the tempo- 
ral domain dj is arbitrary and x represents time. In this 



is the normalized phase- mismatch and j=l,2. 

Physical insight into the properties of Eqs. 111I2[| may 
be obtained from the cascading limit, in which the phase- 
mismatch is large, (3~ x — > 0. Writing E 2 =e 2 exp(i(3z) 
and assuming slow variation of e 2 (x,z) gives the non- 
linear Schrodinger (NLS) equation id z E\ + d\d 2 Ex + 
\E\ \ 2 E\ = 0, in which the local Kerr nonlinearity 
is due to the coupling to the SH field e 2 =E 2 //3. The SH 
is thus slaved to the FW and the widths of the SH and 
FW are fixed. The sign of the mismatch (3 determines 
whether the effective Kerr nonlinearity is focusing or de- 
focusing and thus the cascading limit predicts that both 
bright and dark quadratic solitons can exist. 

However, even for stationary solutions the NLS equa- 
tion is inaccurate, since the term d 2 E 2 is neglected. 
Thus it predicts that in higher dimensions bright soli- 
tons are unstable and will either spread out or collapse 
Q , whereas it is known that stable quadratic solitons ex- 
ist in all dimensions and that collapse cannot occur in 
the x( 2) -system 111 121) 9, 10]. The stabilizing effect of the 
nonlinearity is often described as being due to sat- 
uration of the effective Kerr nonlinearity [jj 0, |n| . We 
show below that the nonlinearity is in fact nonlocal. 

To obtain a more accurate model than that given by 
the cascading limit we assume a slow variation of the SH 
field e 2 (x, z) in the propagation direction only (i.e., only 
d z e 2 is neglected). The SH is still expressed in terms of 
the FW, but now the relation has the form of a convolu- 
tion, leading to the nonlocal equation for the FW 



id z E 1 + d x d 2 x E x + p- l N{El)El = 0, 

/OO 
R(x z)d£, 
-OO 



(3) 
(4) 



with E 2 =j3~ l N exp(i(3z). Equations (l.'ilH clearly show 
that the interaction between the FW and SH is equiv- 
alent to the propagation of a FW in a medium with a 
nonlocal nonlinearity. In the Fourier domain (denoted 
with tilde) the response function R(x) is a Lorentzian 



case d 2 Ej represents pulse dispersion. The parameter (3 R(k) — \/{\ + sa 2 k 2 ), where a~\d 2 / P] 1 ^ 2 represents the 
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degree of nonlocality and s=sign(d2/3)- For s=+l, where 
the x -system 111 121) has a family of bright (for di>0) 
and dark (for e?i<0) soliton solutions 0, R(k) is posi- 
tive definite and localized, giving 

R(x) = (2a)- 1 exp(-\x\/a). (5) 

The cascading limit — > is now seen to correspond 
to the local limit a — ► 0, in which the response function 
becomes a delta function, R(x) — ► 8(x). 

With the nonlocal analogy one can assign simple physi- 
cally intuitive pictures to several important cases. When 
the mismatch \f3\ is large Eqs. (|3I4|) reduce to the so- 
called weakly nonlocal equation with a <C 1 jl3| . Simi- 
larly the nearly phase-matched limit when (3 ss corre- 
sponds to the strongly nonlocal limit with (r»l, when 
Eqs. (i:SI 1H become effectively linear [l4lll5| . 

For s=— 1, R(k) has poles on the real axis and the re- 
sponse function becomes oscillatory in nature with the 
Cauchy principal value R(x)=(2a)~ 1 sin(|x|/<r). In this 
case the propagation of solitons has a close analogy with 
the evolution of a particle in a nonlinear oscillatory po- 
tential. In fact, it is possible to show that the oscillatory 
response function explains the fact that dark and bright 
quadratic solitons radiate linear waves for s=— 1 12] . 

It is important to notice that in contrast to the con- 
ventional nonlocal NLS equation describing, e.g., mate- 
rials with a thermal |l7l or diffusion [l8j type nonlinear- 
ity, liquid crystals [l!j, and photorefractive crystals |20j j. 
the nonlocal response of the x system depends on the 
square of the FW, not its intensity. Thus the phase of the 
FW enters into the picture and one cannot directly trans- 
fer the known nonlocal stability results for plane waves 
[3 EU and solitons [H E3 . 

However, for stationary fields the nonlocal system J3- 
0) represents an exact model for x materials. There- 
fore the known properties of nonlocal solitons in terms 
of profiles directly apply to quadratic solitons. Consider 
stationary solutions to Eqs. 11121) in the form 

E\(x, z) = ai0i(r) exp(i\z), (6) 
E 2 {x, z) = a 2 02(r) exp(i2Az + i[3z), (7) 

where the profile 0y(r) is real, with r = x\\/d\ I 1 / 2 , a\ = 
A 2 1 c^2 / (2di ) | , and a2=A. This scaling reduces the number 
of free parameters to one and transforms Eqs. (11121) into 
the following system 0] 

si<t>'{ - 0i + 0102 = 0, s 2 02 - "02 + 0i/2 = 0, (8) 

where Sj—sigii(Xdj)=±l, a=(2 + f3 / X)\di/ d 2 \, and prime 
denotes differentiation with respect to the argument. The 
properties of solitons described by Eqs. JSJ) are well- 
known • A family of bright (dark) solitons exist for 
S2=si=+1 (s2=— si=+l) and a>0. As discussed above 
we do not consider the combinations s 2 =±s\=— I, for 
which solitons do not exist in the whole a-space. 

The second equation in Eqs. (JSJ) has the formal solu- 
tion 4>2=jN(ip1), with 7=l/(2a) and the nonlocal non- 
linearity N(4>\) given by Eq. For sign(s2a)=+l the 



response function is R(r)=(2a)~ 1 exp(— |t|/ct), with the 
degree of nonlocality a— la]- 1 ' 2 . Inserting the SH into 
the first equation in Eqs. © then gives the exact nonlo- 
cal model for the FW in the system © 

/oo 
fl(r-0$(0d£ = 0, (9) 
-oo 

where 7 is the strength of the nonlocal nonlinearity. Thus 
solitons are equivalent to nonlocal solitons. 
In the weakly nonlocal case cr<d (i.e., |a|^l) the re- 
sponse function R(t) is much narrower than the FW in- 
tensity 2 . Taylor expanding 2 under the integral in 
Eq. we then obtain the weakly nonlocal model 

s 1 9 2 1 -0 1 + 7 (0 2 +a 2 a 2 2 )0 1 = 0, (10) 

where 02=7(1 + a 2 d 2 )<fi 2 . This model has exact bright 
soliton solutions for S2=si=+1 and a>0 ^| 

±T = tanh~ 1 ( ( o) + 2CTtan~ 1 (2a- ( o), (11) 

where p 2 =(a\ — 2 )/(a 2 +4a 2 2 ), a 2 =2/j being the max- 
imum intensity of the FW. Exact dark soliton solutions 
exist for S2=— Si=+1 with nontrivial phase profiles |13| |. 

For |a|<Cl the nonlocality is strong, <73>1, and we can 
expand the response function R(t) in Eq. (0J. For bright 
solitons we then obtain the linear equation for the FW 

Sl 9 2 0! - 0i + 7iVR(<r)0i = 0, (12) 

where 4> 2 =^P\R(t). In this eigenvalue problem the FW 
power P\= 2 (r)c?T plays the role of the eigenvalue, 
and bright solitons correspond to the fundamental mode 
of the waveguide structure created by the exponential 
response function. For s 2 =si=+l and a > Eq. i(T2"|) 
has exact bright soliton solutions in the form of the Bessel 
function of the first kind of order 2a [2lJ 

0i (r) = A 1 J 2s [a 2 ^2P 1 R{T)]. (13) 

For the single-soliton ground-state solution P\ is found as 
the first zero of the derivative, J 2S (-\/ct 3 Pi)=0, which as- 
sures that 0i (0)=0. The amplitude A\ is then found from 
the definition of P u giving Aj^Pi/2 - [Pi/ '(tt 2 ct)] 1/2 . 

In Fig.^we show the full width at half maximum of the 
FW intensity 2 of bright quadratic solitons versus the 
phase- mismatch parameter a. The analytical solutions 
obtained using the nonlocal analogy correctly captures 
the increase of the soliton width with decreasing a. The 
nonlocal model elegantly explains this effect: Because 
of the convolution in the nonlinearity in Eq. repre- 
senting a trapping potential or waveguide structure, this 
potential is always broader than the FW intensity profile 
itself, leading to its weaker confinement and larger width 
when the degree of nonlocality increases. The profiles 
shown in Fig.[I]further illustrates the excellent agreement 
of the numerical results and approximate nonlocal ana- 
lytical solutions in both the weakly (a ^> 1) and strongly 
(a -c 1) nonlocal limit. 
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FIG. 1: Top: Numerically found FWHM (eft) of bright 
quadratic solitons versus a (solid), and the weakly nonlo- 
cal (dashed), strongly nonlocal (dotted), and cascading limit 
(chain-dashed) predictions. Bottom: Numerically found pro- 
files (solid) and strongly nonlocal (left: a=0.01) and weakly 
nonlocal (right: a=10) solutions (dots). S2=si=+1. 

The linear Eq. i|12|) describing the strongly nonlocal 
limit further predicts the existence of multi-hump bright 
solitons. Choosing Pi as the iV'th zero of the deriva- 
tive, i.e., the iV'th root in the equation J 2 ' ff (\/o' 3 Pi)=0, 
gives solitons with an odd number of humps (27V — 1), 
as discussed in However, this does not exhaust all 
soliton solutions supported by the model 1|12|) . There 
also exist antisymmetric solitons with an even number 
of intensity peaks. If the power Pi is found as the TV'th 
zero of the Bessel function itself, and not its derivative, 
i.e., as the TV'th root in the equation J2aW<r 3 Pi)^0 (so 
that 0i(O)=O), then antisymmetric solitons with an even 
number of intensity peaks (27V) exist with the form 

Mr) = s T A 1 J 2(T [a 2 ^2P 1 R(T)l (14) 

where s T =sign(r). When P\, e.g., is fixed by the first zero 
of J 2 a then the solution (|14|) is a two-peak antisymmet- 
ric soliton, which can be interpreted as a bound state of 
two out-of-phasc fundamental solitons. In Fig.[3]we have 
shown the bound state of two out-of-phase fundamental 
solitons predicted by Eq. Ijl4jl . and the corresponding nu- 
merically found solution for a=0.001. We see again that 
the strongly nonlocal model provides an excellent predic- 
tion of this novel bound state quadratic soliton solution. 

In fact, all higher order solitons can be thought of as 
a bound state of a number of individual solitons. For- 
mation of such bound states follows naturally from the 
nonlocal character of the nonlinear interaction. Consider 
two out-of-phase solitons, for which the intensity in the 
overlapping region is always zero. In local Kerr media 




Transverse coordinate t 

FIG. 2: Numerically found bound state of two out-of-phase 
bright solitons for a=0.001 (solid) and the predicted strongly 
nonlocal solution 11411 (dashed). S2=Si=+l. 

the nonlinear change in the refractive index is decreased 
in the overlap region, as compared to the index change 
generated by a single soliton. This leads to a mutual 
repulsion of the solitons. The nonlocality tends to in- 
crease the nonlinear change of the refractive index in the 
overlapping region, and for a sufficiently high degree of 
nonlocality, the index change may even be higher than 
for a soliton in isolation, despite the solitons being out- 
of-phase. This creates an attractive force and leads to 
formation of the bound state. 

In the strongly nonlocal limit the bright soliton is the 
fundamental mode of the waveguide structure R(t) and 
much narrower than the waveguide. In contrast the dark 
soliton is a first mode of the waveguide R(t) at the cut- 
off and as such its width is comparable with that of the 
waveguide. Hence the expansion procedure leading to 
Eq. Ijl2(l is no longer justified. One can use it however, 
but for illustrative purposes only, to show what type of 
dark soliton solutions can be expected in the strongly 
nonlocal system. The linear equation for dark solitons in 
the strongly nonlocal limit will then have the form 

a 1 8^4> 1 -4> 1 +f[Al-Q l R{T)]<h. = 0, (15) 

where Qi = j'^ 00 [A\ — 4>\(r)\dT is the complementary FW 
power and <f)2=l[A\ — Q\R(t)]. For S2=— s%=+l and 
a > Eq. (|15fl has exact dark soliton solutions in the 
form of the zero'th order Bessel function 

Mr) = s T V2^J [a 2 ^2QxR(T)]. (16) 

For the fundamental single-soliton solution Qi is found 
as the first zero Jo[^/o^Qi}=0, which gives Qi = 5.8/ct 3 
and assures that </>i(0)=0. As for bright solitons, choos- 
ing the TV'th root gives novel multihump dark solitons 
with 2N — 1 dips in the intensity profile. However, the 
background amplitude Ai = 2a is fixed and does not sat- 
isfy the self-consistency relation Qi = ^ oa [A\ — </>i(t)]c?t. 

In Fig. |3 we show the full width at half maximum of 
the FW intensity 4>\ of dark quadratic solitons versus 
the mismatch parameter a. The dark solitons have the 
constant background 4>\(±oo) = 2a, ^(ioo) = 1. The 
analytical weakly nonlocal dark soliton solution exist for 
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FIG. 3: Left: Numerically found FWHM(^) of dark solitons 
versus a (solid) and the weakly nonlocal (dashed) and cas- 
cading limit (chain-dashed) predictions. Right: Numerically 
found profiles (dots) and weakly nonlocal solutions (solid) for 
a=10. <j)\ (±oo) = let, 02(±oo) = 1, and S2=— si=+l. 

a > 2 and was taken from Ref. 01 . Unlike bright soli- 
tons, whose width is a monotonic function of a, dark 
solitons are seen to have a clear minimum width of the 
FW at a=a « 3.1. 

Figure [3] confirms that the weakly nonlocal model 
correctly predicts how the soliton width decreases for 
a > ctQ when the mismatch parameter decreases. This, 
as well as the appearance of the minimum in the soli- 
ton width, is again elegantly explained by the nonlocal 
analogy: Because of the convolution in the nonlinearity 
in Eq. @ , representing the trapping potential or waveg- 
uide structure, the contribution from the constant back- 
ground tends to contract this potential. This leads to a 
stronger confinement and thus a smaller width of the soli- 
ton. However, this is only true as long as the amplitude 
of the trapping potential is not affected by nonlocality, as 
in the weakly nonlocal regime. For a high degree of non- 
locality (i.e., smaller value of a) not only the width of the 
trapping potential, but also its amplitude is affected. In 
this regime the nonlocality leads to a drop in the ampli- 



tude of the potential, resulting in a weaker confinement 
and an increase of the soliton width. The profiles shown 
in Fig. |21 further illustrates the excellent agreement of 
the numerical results and approximate nonlocal analyti- 
cal solutions in the weakly nonlocal limit a>l. 

In conclusion, we have used the analogy between para- 
metric interaction in quadratic media and nonlocal Kerr- 
type nonlinearities to provide a physically intuitive the- 
ory for quadratic solitons, which allows for a deeper phys- 
ical insight into their properties and an exact description 
of their profiles. It is well-known that quadratic solitons 
may be characterized by a single solution parameter a, 
which is an effective mismatch parameter depending on 
both the real phase-mismatch and the power. We have 
shown that the nonlocal theory provides a simple and el- 
egant explanation for how the soliton width depends on 
the mismatch a. 

In particular the nonlocal analogy provides simple 
physical models in both the large mismatch limit oC^>\ 
and the near cut-off limit |a|<Cl, corresponding to the 
regimes of weak and strong nonlocality, respectively. 
Our results show that the weakly nonlocal approxima- 
tion gives an accurate description of quadratic solitons 
in a relatively broad range of their existence domain. 
Also, the simple linear physics of the strongly nonlocal 
limit has enabled us to find novel bound states of bright 
quadratic solitons and explain their formation using the 
natural concept of the nonlocality-based attraction be- 
tween out-of-phase constituent solitons. 
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